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Abstract. We study the bosonic Boltzmann-Nordheim kinetic equation, which de- 
scribes the kinetic regime of weakly interacting bosons with s-wave scattering only. 
We consider a spatially homogeneous fluid with an isotropic momentum distribu- 
tion. The issue of the dynamical formation of a Bose-Einstein condensate has been 
studied extensively. We supply here the completed equations of motion for the cou- 
pled system, the energy density distribution of the normal fluid and the density of 
the condensate. With this information the post-nucleation self-similar solution is 
investigated in more detail than before. 
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1 Introduction 



The dynamics of weakly interacting quantum fluids is governed, in approximation, 
by the Boltzmann-Nordheim kinetic equation. There is no a priori restriction on 
either the density or the temperature of the quantum fluid. Since, as explained in 
every textbook on Statistical Mechanics, an ideal Bose fluid undergoes a transition 
from a normal fluid to one with a condensate component, one would expect that 
the kinetic equation retains some information on how the condensate is formed 
dynamically. 

In general, the appropriate kinetic equation is fairly complex and to analyse the 
details of the transition will be a difficult task. There is however one particular 
case which has been investigated in considerable depth. It is assumed that the 
Bose fluid is spatially homogeneous and the momentum distribution is isotropic. 
The Bose particles are massive and interact only through s-wave scattering. The 
distribution function /, as governed by the kinetic equation, depends then only on 
the energy e, e > 0, and on time t, t > 0. Based on theoretical considerations 
and numerical simulations the following scenario has been developed for the kinetics 
of the condensation process. One starts with an initial distribution f{e,0) = f{e) 
which has the density p = dey/ef{e). p is assumed to be supercritical, p > pc- 
Then the solution f{e, t) of the kinetic equation has a piece which concentrates near 
e = and nucleates the condensate at some flnite time t*. More precisely, 

VEf{e,t) = ^f^^^{e,t)+n{t)6{e) (1.1) 

with n(t) = for t < and n{t) > for t > t*. Here fj-^g is a density without any 
delta functions. In the limit t — > oo, n(t) tends to p — pc while the regular piece tends 
to the critical Bose-Einstein distribution (e^^ — 1)^^ with the inverse temperature (3 
determined through the initial energy e = dey/eef{e). Near the solution has 
a self-similar structure, which will be explained below. 

In our contribution we investigate, in more detail than previous studies, the 
precise mechanism of how the condensate is generated and annihilated (on the level 
of the Boltzmann-Nordheim kinetic equation). In particular we obtain additional 
information on the self-similar structure of the solution for t > t^. A numerical check 
of these predictions would help to further elucidate the kinetics of the Bose-Einstein 
transition. 

2 The Boltzmann-Nordheim equation for the en- 
ergy distribution 

We use throughout dimensionless variables in units which minimize the number of 
prefactors. The physically correct dimensions are supplied, e.g., in p]. Instead of 
denoting the energy by e it will be more convenient to use the neutral x, y, z. The 
distribution function is denoted by f{x), / > 0, with volume element ^/xdx, x > 0, 
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since for e = one has d^p = d^lp'^d\p\ and p'^d\p\ = \/2ede in three dimensions. 
To / we associate the density p and energy e as 



/'OO /"OO 

P(/) = / dx^f{x) , e(/) = / dx^fx: 
Jo Jo 



xxf{x). (2.1) 

Physically p is the mass and e the energy per unit volume of the spatially homoge- 
neous Bose fluid, but we stick to the more colloquial expressions. 
fit) is governed by the kinetic equation 

jJ{t)=C,{f{t)). (2.2) 

C4 is the collision operator which describes the scattering of two incoming to two 
outgoing particles, the 4-wave scattering in wave turbulence. For bosons, and under 
the assumptions stated in the introduction, one obtains 

cum = ^ / dydzl{x,y,z){f{x)f{w)f{y)f{z) 

Jv(x) 

-f{x)f{w)f{y)f{z)) . (2.3) 

Here f{x) = 1 + f{x), w = y + z — x, and the domain V{x) = {y, z\y > 0, z > 0,x < 
y + z} for a; > 0. The integral kernel / results from working out the (5-functions 
for energy and momentum conservation of the full 3D collision operator, where x, w 
are the energies of the incoming and [y, z) the energies of the outgoing particles, see 
[H [2] for details. / is defined by 

y/w foTy<x,z<x,0<w, 
Ti \ I \fx ioTx<y,x<z, 

I(x,y,z) = { V r Z Z (2.4) 
^ ^ ^ ^ ioxx<y,z<x, ^ ' 

y/y for ?/ < X , X < 2 . 
If p(/) < 00, e(/) < 00 and if /(t) is bounded, then one can easily work out that 
d 



^-p(/(t)) = / dx^C^U{t)){x) = , (2.5) 







/■oo 

-e(/(t)) = y^ dxV^xC,{f{t)){x)=0. (2.6) 



dt' 

Thus mass and energy are conserved. As we will see, if / which diverges at x = 0, 
the conservation laws may break down. 

On the kinetic level the entropy per unit volume is the one of a non-interacting 
Bose fluid and hence given by 



POO ^ _ 

sif)= dxv/^(/(x)log/(x)-/(x)log/(x)). (2.7) 
^0 
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From (I22D it follows that 

j^sifit)) = aim) (2.8) 

with the entropy production 

ct(/) = / dxdydzl{x, y, z) {A - B) log{A/B) , (2.9) 

J 0<x,y,z<oo,x<'y+z 

A = f{x)f{w)f{y)f{z) , B = f{x)f{w)f{y)f{z) . (2.10) 
Clearly a{f) > 0. 

For a stationary solution cr(/) = 0. Introducing ip = log(///), the condition 
cr(/) = is equivalent to being a coUisional invariant, to say 

ipix) +ip{y + z - x) = ipiy) +ip{z) (2.11) 

for all x,y, z > such that x < y + z. (12.111) admits as only solutions the affine 
functions ip{x) = co+cix. Together with the condition / > this yields as stationary 
solutions of the Boltzmann-Nordheim equation the Bose-Einstein distributions 

ffsAx) = {e^^''^^-l)-\ /3>0,/i<0, (2.12) 

parametrized by the inverse temperature f3 and the chemical potential /x. 
At given (3 the maximal density is 



Pc 



roo poo 

/ dx^/^{e^'' - 1)"^ = , Po= t^xv^(e^ - 1)'^ (2.13) 

Jo Jo 

and correspondingly the maximal energy is 

poo poo 

ec= / dx^/^{e^'' -l)-^x = eo = / dx^/^{e^ - ly^x . (2.14) 

Jo Jo 

Hence there is the critical line 

(e,/eo) = (pc/po)'/' (2.15) 

which divides the (p, e) quadrant into the domain -Dnor of normal fluid and the 
domain Dcon with some fraction of condensate. If, for given initial /, (p(/),e(/)) G 
Dcon, then p(/) > pc and f(t) has no stationary solution to approach in the limit as 
t ^ oo. Based on equilibrium statistical mechanics, it is thus natural to extend the 
class of stationary solutions to 

Vxff3,i,{x) + n^oJix) , neon > , (2.16) 

with the condition that for neon > necessarily p = 0. The density of this distri- 
bution is denoted by p(/3, p) = p(//3,^) for p < and by p(/5, 0, ricon) = Pc + neon for 
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/i = 0. The condensate has zero energy. With this extension there is a one-to-one 
relation <l> between (p, e) and (/?, /x, neon)- 

If in equihbrium there is a (5-function at x = 0, it seems natural to allow for 
such a singular contribution also in the dynamics. The precise formulation will be 
the task of the next section. In fact, Lu [3J proves that, for any / with p(/) < oo, 
e(/) < oo, in the long time limit the solution f{t) converges to (12.161) with the 
parameters (/3, /i, neon) determined through the map $ from the initial data. The 
precise statement requires some preparations and is therefore deferred to Appendix 
C. 

3 Coupled equations for normal fluid and conden- 
sate 

If there is a condensate at time t, the full distribution function is 

Vxftot{x,t) = y/xf{x,t) + n{t)5{x) . (3.1) 

We adopt the convention that / always denotes a continuous function on M_|_ = 
(0, oo) such that p(/) < oo, e(/) < oo and with a possible divergence at a; = 0. We 
postulate that (12. 2p . (12. 3p remain valid when / is substituted by /tot- Inserting /tot 
in (12. 3p yields products of the form ///, 5//, 55/, 555, f f, 5f, and 55. Using the 
continuity of the integral kernel /, it can be shown that the contribution of products 
with more than one 5 vanishes [Hll]. Therefore (12.20 turns into a coupled system of 
equations for / and n, apparently first derived in [5], which reads 

= C,{f) + nCsif) , (3.2) 
-n = —n I dx\/xC^{f){x) . (3.3) 



dt jQ 

The collision operator C3 is familiar from 3-wave interactions. Here it arises since 
one of the collision partners is at rest. Using isotropy and integrating over the 
momentum and energy delta functions, one obtains, 

C,{f){x) = 4= r dy{~f{x)f{x - y)f{y) - f{x)f{x - y)f{y)) 



+ ^ / dy{f{x)f{y - x)f{y) - f{x)f{y - x)f{y)) , (3.4) 

see [1] for details. 

Of course, to properly justify (13.20 . (13. 3p . one has to go back to the microscopic 
model of a fluid of weakly interacting bosons. One imposes an initial quasifree state, 
(■), such that 

{a{k)) = a5{k) , {a{k)a{k')) = 5{k)5{k')a^ , 

{a{kya{k')) = 5{k - k') {\a\^5{k) + /(fc^)) , (3.5) 
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where a{k), a{k)* are the bosonic annihilation and creation operators labeled by mo- 
mentum A; G M^. a G C and |ap is the condensate density. The state (■) is invariant 
under spatial translations. If the BBGKY hierarchy is quasifreely truncated at the 
sixth order, then at least formally one arrives at the coupled system (13. 2p . (13.31) . 

C3 conserves energy, but not mass. By considering the H-Theorem as in (12. 8p . one 
concludes that a{f) = is equivalent to ip = log(///) being a collisional invariant, 
in the sense that 

i^ix) + i^iy) = ^{x + y) (3.6) 

for all x,y >0. Clearly, ip must be linear. Hence the stationary solutions for C3 are 
given by fis^ with /3 > 0, where f3 is determined by the initial energy e(/). 

The coupled system (13. 2p . (13. 3p has an obvious defect. If n(0) = 0, then n{t) = 
for all t, no condensate is nucleated. A second defect lies in the definite sign of the 
term on the right hand side of (13.31) . If / is bounded, then dxy/xC^[f)[x) > 0, as 
will be shown below. Hence, if present initially, the condensate density is monotone 
decreasing. The derivation of (13. 2p . (13. 3p provides no immediate indication of how 
the condensate is generated. Only, since physically the total mass is conserved, 
the mass lost by the normal component has to be identified with the mass of the 
condensate. 

In view of (12.51) valid for bounded /'s a natural guess is that the loss of mass for 
(13. 2p is linked to the divergence of / at x = 0. The relevant properties are stated in 
two propositions, which will be proved in Appendix A. 

Proposition 1. Let f be continuous on M+, p(/) < 00, e(/) < 00, and let 

limx'^^'^fix) =b, b>0. (3.7) 

Then 

poo 

lim / dxy/^Ci{f){x) = -T^b^ . (3.8) 

Here = 3.05 and the defining integral is provided in /IA.13\} . 
Proposition 2. Let f be continuous on ]R_|_, p(/) < 00, e(/) < 00, and let 

limx/(x) = a, a>0. (3.9) 



Then 



roo roo 

lim / dxy/xCs{f){x) = —Tsa"^ + 2 / dxxf{x). (3.10) 
•^^o J 5 Jo 

Here = -2 dxx~^ log(l - x) = tt'^/S. 

Remark: In agreement with the fact that the condensate has zero energy, under the 
condition (13.71) energy is still conserved in the sense that 

roo 

lim/ dxy/xxCj{f){x) = forj = 3,4. (3.11) 
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o 

In f l3.1Up there is a loss of mass as T^a"^ and a gain given by the integral expression. 
We interpret this as the corresponding loss/gain of the condensate. Therefore (13. 2p 
and (13.31) are improved to 

j^m=C,{m)+n{t)Cs{f{t)), (3.12) 

—n{t)=n{t)[T^a{tf-2j dxxf{x,t)). (3.13) 
Note that for ff^^ one has 

POO 

a = P'\ 2 dxxfp^oix) = p-^T^ . (3.14) 

^0 

Thus the Bose-Einstein distributions of (I2.16P are indeed stationary solutions for 

dsn, dsn. 

In ( 13.131) we did not include the term T4^b(tY, since 6 > implies a = oo. Mass is 
transferred between normal fluid and condensate only through nC^. Thus it seems 
that the coupled system (13.120 . (13.130 still has the defect that, in case n(t) = 0, no 
condensate can be nucleated. We will see in the next section, how this objection is 
met. 

Tentatively we assume that a solution f(t),n{t) to (13.12^ . ( 13.13^ has the prop- 
erties: 

(i) fit) > 0, n(t) > are continuous in t, p{f{t)) < oo, e(/(t)) < oo. 

(ii) f{t) : ]R+ ]R_|_ is continuous and limx^Q x f {x , t) = a{t) exists with a{t) < oo. 



4 Scaling theory close to nucleation of condensate 

We note that if n(to) > for some to, then n{t) > for all t > Iq. Thus, if n(0) = 0, 
there must be a first time, t*, such that n{t) = for t < and n{t) > for t > t^. 
Of course, depending on the initial /, one could have = oo. In this section we 
plan to study the solution close t*, the time of nucleation of the condensate. We 
mostly follow [6], but the post-nucleation seems to be novel. 

(i) n{0) = 0, < oo, pre-nucleation of the condensate. 

For the scaling theory the quadratic terms of C4 can be neglected. Physically 
this corresponds to the semiclassical approximation. Since n{t) =0 for t < t^,, we 
have to study 

jtf = c^M (4.1) 

with the semiclassical collision operator 

C4,sc(/)(a:) = ^ / dydzlix, y, z) {f (x) f (y) f (z) + f{w)f{y)f{z) 

V X Jv{x) 

-f{x)f{w)f{y) - f{x)f{w)f{z)) . (4.2) 
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Close to a; = the scaling ansatz is 



f{x,U-t) = -^<P^{-^), 0<t«l. (4.3) 
Inserting in (ji?T|) . (K2\i yields 

r(t)-^-^^r(i)(z/0_(x) + = t-'^+'C,,M.){x) . (4.4) 

We fix the time scale by setting t'^'^'^t = 1. Then, for z/ > 1, 

r{t) = {2{u - l)tY^'^'-'\ (4.5) 

The scaling function becomes then the solution of the nonlinear eigenvalue prob- 
lem ^ 

U(j)^{x) + X — </)_(x) = C4,sc(</>-)(x) . (4.6) 

For r <C X ^ 1 the solution in (14.31) is assumed to be frozen at a definite power 
law. This yields the right boundary condition 

4>-{x) = x'" for a; — > cx) . (4.7) 

For small x the, conventional, left boundary condition is 

0_(x) 00 for X . (4.8) 

In the limit x 0, (14. 6 P then yields, noting that the terms in C4 sc(0+) proportional 
to 00 cancel each other, 

^(t>o= dy dz(f)4y)(f)4z)(l)^{y + z) , (4.9) 
Jo Jo 

which is consistent in the sense that both sides have the same sign. 

p and 00 are free parameters. To have condensation at all, by Proposition 1 
necessarily 

z/>7/6. (4.10) 

On the other hand to have a finite mass in the interval < x < r(t) yields the upper 
bound 

z/<3/2. (4.11) 

Because of slow decay, the scaling function 0_ has infinite mass and energy. It would 
be nice, if (14.61) would determine a unique value of u. But very little is known in 
this direction. 

Numerically one finds that 

= 1.234 (4.12) 
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by fitting the power law of f{x, t^—t) for x > r, see [5l[7j. One also finds /(O, t*— t) = 
T{t)''^(j)o consistent with the right boundary condition fl4.8l) . We refer to [Bj for 
additional information on the eigenvalue equation fl4.6p . 

At the onset of nucleation one has = x~'^, which is outside the space of 

solutions. As we will argue next, for any t > the solution falls back into its proper 
space. 

(a) n(0) = 0, < oo, post-nucleation of the condensate. 

As for t < t^, we can work with the semiclassical approximation, which amounts 

to 

^/ = C4,sc(/)+^C3,sc(/), (4.13) 

^n = nT3a'^. (4.14) 
The 3-wave semiclassical collision operator reads 

C3,sc(/)(x) = ^ r dy{f{x - y)f{y) - f{x)f{x ~ y) - f{x)f{y)} 

+ ^ / dy{f{x)f{y) + f{y)f{y - x) - f{x)f{y - x)} . (4.15) 



'X _ ^ 

We impose the scaling ansatz for x close to as 

/(a;,t, + t) = -ip0+(-|y), 0<t«l. (4.16) 

The boundary condition at infinity is 

(f)+{x) = x^'^ for X ^ oo . (4.17) 
If / denotes the continuation of the scaling solution from (i), then one should 

set 

u = iy, (4.18) 

since for a very short time span and away from x = the distribution function is 
frozen. On the other hand a general u is of interest, since one can set = and 
impose by hand some power law for the intial distribution, i.e. one sets 

f(x) = x-^h{x) for t = (4.19) 

with a cutoff function h so to have bounded p(/), e(/). Here necessarily 

l<z><3/2, (4.20) 

the lower bound being required for nucleation, the upper bound for p(/) < oo. 
The condensate density is assumed to scale as 

n{t) = {t/toV. (4.21) 
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Inserting in (14.140 yields 

7t-i = T,a{tf . (4.22) 
To ensure a > 0, we have to impose the left boundary condition 

(f)Jx) = — for X ^ , ao > . (4.23) 

X 



Then 

ait) = 

x—>Q ' ' ^ T 

which implies 



a{t) = lim XT-y+{-) = T-^+^ao , (4.24) 



r(t)-^+i = -(-^)i/Vi/^ (4.25) 

Next we insert the scaling ansatz fl4.16p . f l4.2ip in (14.131) . Multiplying both sides 
by T^'^~'^ one arrives at 

- r2--3f(z70+ + x^;) = C4,sc(0+) + (tAo)^r^"'/'C3,sc(0+) . (4.26) 

The prefactor on the left turns out to be independent of t. Balancing the prefactor 
of C3,sc yields 

7 > provided 1 < u < 3/2 consistent with (I4.20p . We conclude that the scaling 
function satisfies 

- K,{U^+ + X^'^) = C4,sc(0+) + i^2C3,sc(0+) (4.28) 

with the boundary conditions (14.171) and (14.231) and the two positive constants 

K, = ain{i-ur\ K2 = {-^^y. (4.29) 

Co and to are free parameters. 

There is one important consistency check. In (I4.28P we take the limit x — 0. 
The left hand side behaves as —Kiao(u — l)/x. For C^^sc we use, see Appendix A, 



/•OO /"OO ^ 

lim xC4,sc (</>+) (a;) = (ao)^r4,i , Ti^i = dy dz{{y + z - l)yz) . 

Ji Ji 



(4.30) 



Numerically r4i = 1.645. Clearly, one cannot balance the negative —KiaQiv — 1). 
Therefore, C^^sc must come into play. For small arguments the leading terms in the 
integrand of C^^sc cancel and, to achieve the 1/x singularity, we have to assume the 
subleading behavior 



X 



-1 = aox + aix^/^ for X . (4.31) 
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Then, see Appendix A, 



lim xC' 



■3,sc 



(ao)^air3,i , 



r 



3,1 



f 1 

/ du— -(2 

Jo u{l-u) 



3/2_(i_^)3/2)^ 



(4.32) 



Numerically Fs^i 



5.56. We arrive at 



Kiao{u - 1) 



aor4^i — K2aQaiTs,i , 



(4.33) 



which determines Oi. Thus 0+ has a subleading behavior as — aix~^/^ for x — >■ 0. 
By the same procedure one could, in principle, compute the next order corrections 

to 0+. 

At first sight the nonlinear eigenvalue problems (14.61) and (I4.28P look rather 
similar, except for a further collision operator and the change of sign on the left. 
On physical grounds we expect however a very different solution behavior. (14. 6 p 
should have a single eigenvalue, or possibly a discrete set of eigenvalues [6j, while 
(I4.28P should have for any u with 1 < u < 3/2 a. solution satisfying the boundary 
conditions (HTT]) and (H:23|) . 

Physically the most interesting, and accessible, prediction is the exponent 7 of 
(I4.27P which governs the initial increase in the condensate density. If in (14.191) . one 
sets 7/6 < z/ < 3/2, then, according to Proposition 1, C4 provides a mechanism for 
nucleation. Immediately, i.e. for any t > 0, the divergence at x = drops to a{t)/x, 
a{t) from (I4.22p . For u = u = 1.234 numerical simulations are available [7] and yield 
7 = 0.571 in good agreement with the prediction (14.271) . At u = 7/6, (14.271) results 
in 7 = 1. Aspects of this case have been established mathematically by Escobedo et 
al. and we state their results in Appendix B. Also the numerical solutions in [H [5] 
clearly show a linear initial increase in the condensate density and the rapid switch 
from the initial divergence x~'^^^ to the slower x~^. For 1 <1j < 7/6 and n(0) = 0, 
no mechanism for nucleation at time t = is available. (Of course, at some later 
time the scenario as described under (i) may set in). Thus to verify the prediction 
(I4.27P one would have to start with some small n(0), so to set C3 in action. 

5 Conclusions 

In the "naive" picture of the condensation process the distribution function / de- 
velops a 5-like concentration close to the origin, either at some finite time or in the 
long time limit. The Boltzmann-Nordheim kinetic equation tells a different story, 
however: If there is sufficient mass accumulated near the origin, the solution de- 
velops explosively a 1/x singularity at (which is integrable for y/xdx). Once the 
singularity is formed, mass of the normal fluid can be channeled into the condensate. 
On the other hand the condensate is annihilated at a rate 2 dxxf{x, t). For long 
times both processes balance so to approach the equilibrium condensate density. 
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A Appendix: Proofs of Propositions 1, 2 and Eqs. 

(USUI), 032). 



We start with the shghtly easier proof of Proposition 2. 
Proof of Proposition 2: We have to study the hmit 5 ^ of 

/•oo r-x 

2 dx dy{{fix - y)f{y) - f{x)f{x - y) - f{x)f{y)} - fix)) 
Js Jo 

POO poo 

+4 y dxj dy{{f{x)f{y) + f{y - x)f{y) - f{x)f{y - x)} + f{y)) 
= Ar{5) + M5) + A:,{5) , (A.l) 

where Ai{5) is the integral with the first curly bracket, ^2(5) the one with the second 
curly bracket, and A^{5) is the sum of the two integrals linear in /. 

We consider Ai[5). Each summand is linearly transformed such that the in- 
tegrand is of product form, f{x)f{y). For /(x — y)f{y) the new domain of in- 
tegration is {x,y > 0\6 < X + y}, for f{x)f{x — y) the new domain of integra- 
tion is {x^y > 0|x < y, S < y}, and for f{x)f{y) the domain of integration is 
{x^y > 0\6 < X , y < x}. Altogether 

A^i6) = 2 f dx f dyf{x)f{y). (A.2) 

Jo Js-x 

By continuity one has the bounds a-{6)x~^ < f{x) < a+{6)x—l valid for < a; < 5 
with lims^o a±{S) = a. Inserting in ( ]A.2I) yields 

Tsa46f<A,{6)<Tsa+{6)\ (A.3) 

hence 

\im Ai{5) = T3a\ (A.4) 

Next we study A2{S) by the same technique. For f{x)f{y) the domain of in- 
tegration is {x,y > 0\S < X < oo,x < y}, for f{y — x)f{y) the new domain of 
integration is {x,y > 0|0 < x < oo,x + 6 < y}, and for f{x)f{y — x) the new 
domain of integration is {x,y > 0\6 < x < oo}. Altogether, using the symmetry of 
the integrand, one obtains 



oo 



A,{5) = -A dx dyf{x)f{y). (A.5) 

Js Jx-S 

We split the x integration into 6 < x < 6i and 6i < x < oo, with 6 < 6i. For the 
second interval 

/"OO pX 

hm / dx / dyf{x)f{y) = (A.6) 

J&x Jx-S 
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by dominated convergence, since dxf{x) < oo and f^^_^dyf{y) ^ for 5 ^ 
with di < X. For the interval 6 < x < 6i we use the bounds a^{di)x^^ < f{x) < 
a+((5i)x^^ vahd for < a; < 5i with lims-^^o a±{6i) = a. Letting 5 — > and 
subsequently (5i — > 0, we conclude 

\imA2(S) = -2Tsa\ (A. 7) 

5-»0 



Finally for the linear terms we use that 

/■oo 

/ dxf{x)x < oo . (A.8) 
Jo 

Therefore 



poo 

lim (Ai{5) + A2{5) + AsiS)) = -Fga^ + 2 / dxxf{x) . (A.9) 

Jo 



Proof of Proposition 1: We have to study the limit 5 ^ of 



□ 



dx / dydzl{x, y, z) {{f {x) f {y) f {z) + f{w)f{y)f{z) - f{x)f{w)f{y) 

Jv{x) 

-f(x)f(w)f(z)} + f(y)f(z) - f(x)f(w)) 

^ 3,(5) + B^iS) , (A.IO) 

where Bi{S) is integral with the curly brackets and 82(6) is the one quadratic in /. 
We consider first Bi{5). For each summand we make a linear change of variables 

to the product form f (x) f (y) f (z) . For f{x)f{y)f{z) the domain integration is 
{x.y.z > 0\d < x.x < y + z}. For f{w)f{y)f{z) the new domain of integration is 
{x, y,z > 0\x + S < y + z}. Using the symmetry of the integrand, for f{x)f{w)f{y) 
the new domain of integration is {x,y,z > 0\x < y + z, z > S} andfor f{x)f{w)f{z) 
the new one is {x,y,z > 0\x < y + z,y > 5}. Under these transformations the 
integral kernel / is not altered. Adding the four terms yields 

r><5 px+S fX—y+S 



po px+0 px—y+d 

B,{5) = -2 / dx / dy / dz^zf{x)f{y)f{z) 
Jo Js Jo 

rS pS pS 

+ dx dy dzy/xf{x)f{y)f{z) 

Jo Jx Jx-y+5 



f 

Jo 



no poo poo 

i dx dy dz^/xf{x)f{y)f{z) 
Js Js 

00 px px—y+S 

dx I dy I dz^y + z- xf{x)f{y)f{z) 

pco pS poo 

/ dx dy dz^f{x)f{y)f{z). (A.ll) 

JS Jo Jx-v+S 
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As in the proof of Proposition 2, one sphts the domain of the x integration into < 
X < 5i, 5i < X < oo with 5 < 5i. By continuity one has the bounds < 
f{x) < yahd for < X < 5i with hm^^^o &±(5i) = b. Since p(/) + e(/) < 

oo, we conclude that in the limit 6 —>■ the integral over the interval 6i < x < oo 
vanishes. Inserting the bounds, the integral over < x < 6i becomes 

T4b^{6if < \imBi{6) < ^46+((5l)^ (A.12) 

(5— >0 

where r4 is independent of 6i and given by 

pi pl+x px—y+l 

-T4 = -2 dx dy dz^/^{xyz)-^/^ 
Jo Jl Jo 

+ dx dy dzy/x{xyz)^'^^^ 

Jo Jx Jx-v+1 



' x—y+1 

'1 POO POD 



PI POO POO 

— dx dy dzy/x{xyz)~'^^^ 
Jo Jl Jl 



00 px px—y+1 

dx dy dz\/y + z — x{xyz)~'^^^ 

1 Jl J rciax(x~y,l) 

POO pi POO 

+2 dx dy dz^{xyzY'^l^ . (A.13) 

Jl Jo Jx-y+l 

The result follows by taking 61 0. 

Finally we discuss B2{S). Under our assumptions each one of the integrals is 
finite separately. Therefore 

POO P 

limi?2(5)=/ dx dydzIix,y,z){fiy)f{z)-fix)f{y + z-x))=0,iA.U) 

■5^0 Jo Jv^ 

by a linear change of coordinates in the second summand. □ 

r4 cannot be computed explicitly. Numerically one finds r4 = 3.05. 
We supply the proofs of fOO]) and (gSD- 

Lemma 3. Let f be continuously differentiable on R+, / > 0, dxf{x) < 00, 
and 

lima;/(x) = ao > . (A. 15) 



x^O 

Then 



POO POO 

lima;C4,sc(/)(a;) = / dy dzi{y + z - l)yz)-K (A.16) 

Jl Jl 

Proof: The domains in the definition of V{x) are labelled by I, II, III, IV in the 
order of (12. 4p . Each domain will be discussed separately. 
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ad I: All arguments are small and one can use that f{x) = a^/x for small x. Then 
one has 

-1 px px 

xf{x)^ dy dzV^{f{y)f{z)-f{w)f{y)-f{w)f{z)) 

"I PX P X 

+— dy dz^f{w)f{y)f{z) 

10 Jx-y 



X . 

3 



(flo) 



j dy j dz[yz^Jy + z - l) ^ 

Jo Jl-y 



+ I dy 1^ dz{yzy/y + z-l)-'^ = (A.17) 



valid for small x, where we use that in the last line the integral is well defined. 
ad II: We first note that 

/•OO /"OO 

xfix) dy dz{f{y)f{z)-f{y + z-x)f{y)-f{y + z-x)f{z))=0,{A.18) 



since each term of the integrand is integrable. For the remaining integral we choose 
a fixed 6 > and split the y,z integration into [x,S] and [5, cx)). The integral 
over [S, oo) x [6, oo) is proportional to x and the off-diagonal part is proportional to 
xlogx. Thus it remains to study the domain [x, S] and [x, 6]. Choosing S small and 
approximating /(x) by ao/x, one obtains 

(ao)^a^ dy dz{{y + z — x)yz^ = {a^fx I dy I dz(^{y + z ~ l)yz) 



(A.19) 

which tends to the expression in flA.161) as x — > 0. 

ad III (and ad IV by symmetry): We fix 5 > such that {xf{x) — < c{6) with 
c(5) — > as (5 — > 0. We split the integral over y in [0,6] and [5, oo). For the first 
domain the argument in ad I shows that the limit x ^ vanishes. The remaining 
integral is 



X I dy I dz^z{f{x)f{z)U{y)-f{y + z~x)) 

+f{w)f{y)f{z) - f{x)f{w)f{y)) (A.20) 

POO PX 

X dy dzy/z{aoX~'^f{z)f'{y){x - z) + fiyfa^z"^ - aox'^f^yf) , 



which is of order x and vanishes as x ^ 0. 

Lemma 4. Let f be once continuously differentiable on M.^, f > 0, dxf{x) < oo 
and there exists constants e,6, e > 0, 6 > 0, qq, bo, cq, such that 

\f{x)-^ - aox - box^/^l < cox^/^+' (A.21) 
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for X G [0, 6]. Then 

limxC3,sc(/)(a:) = {ao)\ f du (2~Au-'/^){l-u'/'-{l-uf"). (A.22) 

x^O Jq u[l - u) 

Proof: We study both terms of Cs^sc separately. 

Term 1. Since x is small we can use the approximation (1A.21|) to obtain 



x^ r dyfix - y)f{x)f{y) {f{x)-' - f{y)-^ - f{x - y)-') 
Jo 

- {ao)%x^ r dy{{x - y)yxr\x''^ - y'" - {x - yf'^) 
yx Jo 

= 2{ao)\ f dy{{l - y)y)-\l - y'/' - (1 - yf") . (A.23) 
Jo 

Term 2. As before we split the y integration into [x, 5] and [5, cxd). For the first 
integral we follow (]A.23|) and arrive at 

x^ f dyf{x)fiy)f{y - x) {f{y - x)-' + f{x)-' - fiy)-') 

= 4{ao)\ [ ^dy{y{y-l))-\{y-lf/' + l-y'/'), (A.24) 



which tends to the expression in (]A.22p as x — 0. 



For the second integral one has 

/•oo poo noo 

^Vx dyfiy -x)f{y)+Af{x)y/^{ dyf{y)- dyf{y-x)) 
Js Js Js 

POO 

= 4y^y^ dyf{yf + Afix)V^xf{5), (A.25) 
which is proportional to y/x and thus subleading. 

B Appendix: Singular solutions for the Boltz- 
mann-Nordheim equation 

We state in more detail the recent results of Escobedo, Mischler, and Velasquez 
[8]. They consider the Boltzmann-Nordheim equation (12. 2p with initial / which 
diverg ' for X — 0. More precisely they assume that the initial / is once 

continuously differentiable on (0, oo) and satisfies 

\f{x) - Ax-'^/^\ < Bx-^'^/^-^\ 0<x<l, (B.l) 
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\f'{x) + lAx-^^/^\ < Bx~^^^'^-^^ , < X < 1 , (B.2) 

f{x)<Be-^\ l<x (B.3) 

for some positive constants A, B, D, and 6. Note that < oo, e(/) < oo and 
that there is no condition on their size. In particular p{f) could be much smaller 
than pc- 

Theorem (EMV). For any f satisfying ilB. 1\) to W. 3\) . there exists a unique solu- 
tion, f{x,t), to h2.2^) . continuously differentiable in t /or < t < oo and continuous 
in X for < x < oo, as well as a function b{t), satisfying 

< /(x,t) < Lx-^/V^% x>0, tG(0,T), (B.4) 



-6(t)a;"^/*^| < La;-(^/^-^/2) , < a; < 1 , tG(0,T), (B.5) 
0<b{t)<L, tG(0,T), (B.6) 
for some positive constant L and for some T > 0, T depending on A, B, and 6. 
Together with our Proposition 1 it then follows that, for < t < T, 



e(/(t)) = e(/) , (B.7) 
p{f{t))=p{f)-T, fdsb{s)\ (B.8) 







n{t) = T4 / dsb{sy. (B.9) 
Jo 

Thus the condensate density increases linearly in t for small t. 

It would be of interest to extend such a result to the coupled system. As argued, 
we expect that C3 dominates C4 and therefore in fIB.SP the divergence at x = should 



drop from to x~^, as also reported in [IJ for numerical solutions. 

Currently, the restriction to bounded T is needed for the proof. Whether there 
is really such a restriction can only be speculated. One could argue that by the 
choice of initial conditions one has opened already a channel for the formation of 
condensate. If the system uses only that channel, then 

/■oo 

14 / dtbitf < pif) (B.IO) 

and necessarily b(t) — » as t — 00. Once f{t) is subcritical, the mass transfer could 
stop and f(t) should converge for large t to the corresponding unique equilibrium 
distribution. 
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C Appendix: Long time behavior 

Lu [3] considers Equation (12.21) in the weak form, 

-j^ ^{x)f{x,t)y/^dx = J ip{x)C,{f{t)){x)V^dx (C.l) 

with test functions ip G C^(]R+) and proves that (1C.1I) remains meaningful even if 
■\/xf{x, t)dx is substituted a positive measure on IR+. We denote such a measure by 
f{dx,t). The possible roughness of f{dx,t) is balanced by rewriting the right hand 
side of (IC.ip in such a way that the colhsional difference ip{x) +ip{y + z — x) —(p{y) — 
ip{z) appears. For t = we assume finite mass and energy, i.e. f{dx) < oo and 
J^xf{dx) < oo. Then ( 1C.ll) has a solution f{dx,t), which conserves mass and 



energy. It is not known whether the solution constructed by Lu is identical to the 
solution of the coupled system ( 13.121) . ( IS.lSp . The stationary measures of ( 1C.1I) are 
necessarily of the form 

ftot,p,f,{dx) = ff3,f,{x)y/xdx + ncoJ{x)dx (C.2) 

with either /i < and condensate density neon = or /i = and ricon > 0. 

Let us assume that the initial data are given by a density, i.e. f{dx) = fa_c{x)^/xdx, 
such that 



P = P{f) 



POO POO 

/ f^c{,x)y/xdx < oo , e = e(/) = / xfi^{x)^/xdx<oo. (C.3) 
Jo Jo 



In general, the solution to (IC.ip is then a measure, f{dx,t). It can be uniquely 
decomposed into an absolutely continuous and a singular part, 

f{dx, t) = U{x, t)^dx + /s(rfx, t) . (C.4) 

Now let fp^^ be the equilibrium distribution corresponding to p(/), e(/), where /i = 
in case neon = P — Pc > 0. Lu [3J proves that, for p < Pc, 

/»oo /»oo 

lim / f^{dx,t) = 0, lim / \fa.c{x,t) - f|3^f,{x)\^/xdx = . (C.5) 
On the other hand for p > Pc one has 

POO 

lim / \fac{x,t) - ff^^o{x)\x^/xdx = (C.6) 



and 



t~*oo 



lim / f{dx, t) = p- Pc, (C.7) 



' {x\0<x^<r{t)} 

where r{t) is the integral in (1C.6I) . Thus for large times the solution has a 5-peak 
like concentration in the interval [0,r(t)^/^] with weight n^on- 
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